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Abstract. Light radioactive nuclei play an important role in many astrophysical environments. 
Due to very low cross sections of some neutron and proton capture reactions by these radioactive 
nuclei at energies of astrophysical interest, direct laboratory measurements are very difficult. For 
radioactive nuclei such as ^Li and ^B, the direct measurement of neutron capture reactions is 
impossible. Indirect methods have been applied to overcome these difficulties. In this work we will 
report on the results and discussion of phenomenological potential models used to determine some 
proton and neutron capture reactions. As a test we show the results for the ^^0(p,7)'^Fi;s(5/2+) and 
'^0(p,7)^^Ff,x(l/2^) capture reactions. We also computed the nucleosynthesis cross sections for the 
^Li(n,7)^Ligi, ^L,\{n,yflAgs and ^B(p,7)^Qf capture reactions. 
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INTRODUCTION 



In the investigation of many astrophysical entities such as primordial universe, main 
path stellar evolution, novae, super-novae explosion. X-ray bursts etc. , the important in- 
put parameter in the models is the cross sections of the capture reactions as a function of 
energy. However, since in many of these environments the temperature is not very high, 
the cross sections have to be obtained at the "Gamow Peak" energies which are very 
low, in the range of few tens to at most hundreds of keVs. Whereas some of the cross 
sections originate from laboratory measurements, the majority are based on extrapola- 
tions from the higher to lower energies or are obtained by pure theoretical models with, 
sometimes, no firm experimental basis. Some reactions are practically impossible to be 
directly measured. For instance, for light radioactive nuclei nucleosynthesis, the neutron 
capture reaction is impossible because the combination of target-i-beam is not possible. A 
typical case is the ^Li(n,7)^Li capture reaction where direct measurement is not possible 
because no ^Li or neutron target exist. For proton capture reactions involving radioactive 
nuclei such as or ^Be, the direct measurement is very difficult due to the low cross 
sections and limited beam intensities of these elements. For such cases, indirect methods 
to determine the cross sections were developed and are quite usually adopted. Extrap- 
olation from higher energy data to lower regime is also not straight forward. Careful 
and accurate account of physically relevant information has to be considered in the de- 
scription of the reaction before the extrapolation is performed, where to describe capture 
reactions not only information on the structure of the nuclei, but also a clear understand- 



ing of the reaction mechanism is required. 

Many indirect methods have been developed to overcame the difficulties to obtain 
low energies cross sections. Among these methods we have the Coulomb dissociation 
method, which corresponds to the inverse temporal reaction of the capture [1], the trojan 
horse method [2], the reduced- width or ANC (Asymptotic Normalization Coefficient) 
method [3] and potential model, where the latter two use transfer reactions as a way 
to get information on the non-resonant part of the capture reaction process and will be 
discussed in more detail in the next section. These indirect methods are very suitable to 
be used in association with low-energy radioactive nuclear beams. 

LIGHT RADIOACTIVE NUCLEI AND NUCLEOSYNTHESIS. 



To overcame the A=8 gap and synthetize heavier elements the key nuclei are the light 

radioactive elements ^Li(178 ms) and *^B(770 ms). These elements are important for 
instance in the inhomogeneous big-bang nucleosynthesis (IBBN) [4]. Because of the 
assumption of homogeneous baryon density distribution and also the instability of ^Be, 
the nuclear reaction flows stop at A=8 in SBBN (Standard Big Bang Nucleosynthesis). 
In IBBN the baryon density distribution is assumed to be inhomogeneous due to several 
cosmological processes before an onset of primordial nucleosynthesis. A difference in 
diffusion between protons and neutrons makes high density zones proton-rich and low 
density zones neutron-rich. In such an environment, heavy elements can be produced 
via the formation proton or neutron rich unstable nuclei. Also it is believed that in 
supermassive stars with low metalicity, the 3a capture process cannot form enough 
^^C to initiate the rapid hydrogen burning during the explosion. However, in such high 
temperature and densities the hot pp-chmn my bypass the slower 3 a capture [5] and in 
this case ^B plays an important role. 

The following nuclear reaction chains were found to be very important in nucleosyn- 
thesis processes discussed above to jump the A=8 gap and synthetise heavier elements: 

chain (1): ^Li(n,7)Ma,n)iiB(n,7)i2B(e-v)i2c. 

chain (2): ^Li(n,7)^Li(n,7)'^Li(e-v)'^Be(n,7)i0Be(e~v)i^B(n,7)iiB(n,7)i2B(e-y)i2(-._ 
chain (3): ^Be(p,7)^B(a,p)iiC(e+v)iiB(p,7)i2c. 
chain (4): '^Be(p,YfB(a,py^C(p,Yy^N(l5y^C. 
chain (5): '^Be(p,rfB(p,YfC(a,p)^^N(l5)^^C. 

The first two reaction chains play the central role in the production of neutron-rich 
nuclei, and the chains (3), (4) and (5) of proton-rich nuclei. The ^Li(n,7)^Li is also an 
important reaction in the early stage of the ignition of supemovae, where two neutrons 
capture reactions,^Li(2n,7)^^Li and ^He(2n,7)*^He(/3)^Li, may become important. The 
^Be(p,7)^B is related to the solar neutrino problem. Precise predictions of the produc- 
tion rate of ^B solar neutrinos are important for testing solar models, and for limiting 
the allowed neutrino mixing parameters. This, however, it the most uncertain reaction 
leading to ^B formation in the Sun [6]. The ^B(p,7)^C, capture reaction is also impor- 
tant in the novae environment where temperatures are several times larger than 10^ K, 



corresponding to Gamow window energies around £■ = 50 — 300 keV [5, 7]. 



5-FACTOR AT LOW ENERGIES 



Cross sections for many of the neutron and proton capture reactions by light radioac- 
tive nuclei are still poorly known or remain unmeasured at the required low energy re- 
gion. For systems which cannot be measured directly at low energies some procedures or 
methods are adopted to obtain cross sections in this region. The most obvious way to ob- 
tain the low energy cross section is the extrapolation to low energy through polynomial 
parametrization of the available high energy region data. Actually these extrapolations 
are performed to the S-factor where the strong energy dependence of charged particle 
capture reactions due to the coulomb barrier penetration are removed. Thus, the cross 
sections are conventionally expressed in terms of the ^-factor which is defined as: 

S{E)=Eo{E)exp[27i:ri{E)] with ri{E) ^ZaZbe^/hv, (1) 

where v is the initial relative velocity between the nucleus and proton. 

However, the S-factor does not remove completely the energy dependence, the struc- 
ture of the final bound state, resonances and attenuation of the barrier by the nuclear 
mean field can give an extra energy dependency. 

The problem is that the extrapolation is not usually a straight forward procedure and 
some physics has to be taken into account. For example, for some systems there is an 
upturn in the S-factor at very low energy due to a pole at threshold when the photon 
energy of aA{x,Ey)B capture reaction vanishes. This is a quite common feature and a 
straight line extrapolation of S{E) from the high energy data region to S(0) would not 
work. According to Jennings et al. [8], the presence of the pole suggests the 5-factor 
may be parametrized as Laurent series S = dix E~^ -\-d[i-\-d\ x Ej, and they propose 
the following expression to parametrize the S-factor: 

|7^ = Tr^+^ + ^x^ or S{E) = -^+B + CxE (2) 
S(0) Eb-\-E Eb + E 

A good test for this assumption is the investigation of the low energy region of the 
^^0(p,7)^^F capture reaction. This reaction is an important reaction in the CNO cycle of 
our Sun but also in evolutionary phase referred as the asymptotic giant branch (AGB) 
of some massive stars [9]. A very good and precise set of data has been obtained for 
this reaction for both ground-state (5/2+) and first excited state (1/2+) transitions of ^^F, 
as can be seen in Figure 1. In the figure the upturn for the reaction leading to the 1 /2+ 
excited state is very clear. 

Let's consider only the data from Rolfs [11] of the transition to the ^^Fgs(5/2+) 
ground-state. By fitting the data with the expression with pole effect and by a straight 
line we get the following parameters, respectively: 

a) S{E) = -0.4359 + 0.4584 x E, giving S(0)=1.326 and X^=0.194 

b) S{E) = 0.66189 + 0.1 x E, which gives S(0)=0.66 mdx^^0305 




As we can see a very different value for 5(0) is obtained depending on the extrap- 
olation performed. The solution for this problem is including more physics in the ex- 
trapolation to better determine the 5- factor at low energies. This can be better done with 
another ''fitting" method called R-matrix [13]. The idea of the R-Matrix calculations 
is to parameterize the cross sections with a small number of parameters and then used 
to extrapolate the cross section down to astrophysical energies. Since it deals with res- 
onances in the continuum, it can be considered, to some extent, a generalization of the 
Breit-Wigner formalism. The R-matrix framework assumes that the space is divided into 
internal region (with radius match a), where nuclear forces are important, and the exter- 
nal region, where the interaction between the nuclei is governed by the Coulomb force 
only. For more details of this method we recommend the work of R Descouvemont in 
ref. [13]. This method, however, has some limitations. It doesn't work for all systems, 
the radius match a is not clearly determined (usually it is arbitrary chosen) and it is valid 
only for low level density composite nuclei. Moreover, it has adjusted parameters which 
not always has a clear physical meaning. 

Other possibility is to use what we call ''non — Fitting" methods. These other methods 
are related to theoretical calculations and they determine the cross sections from wave 
functions with, in principle, no data is required. Among these methods we have the 
microscopic cluster models and Microscopic "aZ? initio" models. In the first, a many- 
body Hamiltonian with embedded cluster are involved in a more complex calculations. 
Internal structure has to be taken into account and the resonating group method can be 
applied. In the second, even more complicate calculations are involved since continuum 
states, which are very difficult to take into account, have to be considered. Applications 
of these two calculations to the cases of ^Li(n,7)^Li and ^Be(p,7)^B capture reactions 
can be seen in ref. [14]. Another method also considered as "non — Fitting" method is the 
Potential model, where ANC (Asymptotic Normalization Coefficient) can be considered 



as an approximation. This method is quite simple to solve numerically since no structure 
of the colliding nuclei has to be taken into account. 



THE POTENTIAL MODEL 



In the potential model, the direct radiative capture (DRC) of an s- and/or J-wave 
nucleon (proton or neutron) by a nucleus A, proceeding via El transition and leaving the 
compound nucleus B in its ground state, is given by: 

ofi5(n,7) = ^kl\<^,cat\0^'\^bound > |', 0) 

where kj = Ej/hc is the wave number corresponding to a y-ray energy Ej, O^^ stands 
for the electric dipole operator, the initial-state wave function ^scat is the incoming 
nucleon wave function in the nucleon-nucleus potential and ^tomd is wave function 
which describes the single-particle bound state. 

Although the potential model works well for many nuclear reactions of interest in as- 
trophysics, it is often necessary to pursue a more microscopic approach to reproduce ex- 
perimental data. In a microscopic approach, instead of the single-particle wave functions 
one often makes use of overlap integrals, Ihoundi^), and a many -body wave function for 
the relative motion, %caj (r). The effect of many -body will eventually disappear at large 
distances between the nucleon and the nucleus. One thus expects that the overlap func- 
tion asymptotically matches the solution of the Schrodinger equation, with V = Vcoui 
(Pure Coulomb) for protons and V = for neutrons. This assumption may be true only 
for a very peripheral capture reaction. This approximation, when r — > oo, is called ANC 
(Asymptotic Normalization Coefficient) and: 

houndir) = ANC X ^-^-^z^+i/^^^^^) (for protons) (4) 

r 

where the binding energy of the A+ x system is related to k by means of Eb = 
h^K^/lninx, Wp^q is the Whittaker function and ANC is the asymptotic normalization 
coefficient. The ANC can be obtained from peripheral transfer reactions whose ampli- 
tudes contain the same overlap function as the amplitude of the corresponding capture 
reaction of interest [3]. 

In terms of overlap integral the direct capture cross sections are obtained from the 
calculation of 

<4 °^ I < \\r^YL\ l^scatir) > |". (5) 

However, it has been shown that 5-wave neutron capture, even at rather low energies, 
is not peripheral [15, 16] and so it is necessary to calculate the wave function of the 
incoming neutron or proton and the wave function for the bound system. Thus, in the 
potential model, it is necessary to calculate the overlap function also taking into account 
the internal part of the nuclear potential, lbound=SF^/-^ x ^{r) hound- Here, SF^/-^ is the 



spectroscopic amplitude and ^{r) hound is the wave-function which describes the bound 
state. 

In the potential model, the continuum wave function, ^scat, has to be calculated with 
a potential which includes also the nuclear interaction. Thus, the essential ingredients 
in the potential model are the potentials used to generate the wave functions ^scat and 
hound-, and the normalization for the latter which is given by its spectroscopic factor. 
Here we use a Woods-Saxon (WS) parameterization to build up the potentials Vo{r) 
and Vsoir), where for the latter we consider the derivative of the WS form factor. The 
parameters Vq, Vsq, Rq, ao, Rso, and 1350 are chosen to reproduce the ground state energy 
Eb (or the energy of an excited state). For this purpose, we define typical values (Table I) 
for V50. ^0. ^0-, Rso, and vary only the depth of the central potential, Vq. A different set of 
potential depths might be used for continuum states. To calculate the non-resonant part 
of these capture reactions in the framework of the potential model we used the computer 
code RADCAP developed by Bertulani [17]. 

We have tested this potential model by determining the ^^0(p,7)^^Fgi, and 
^^0(p,7)^^Fij^ ^-factor as a function of energy. In Table-I we list all the parame- 
ters of the potentials used to generate the incoming and bound wave functions. All 
the potentials were assumed to be a Woods-Saxon shape with geometric parameters 
ro = 1.25 and a — 0.65 fm. The Jb — 5/2^ ground state {Jb =1/2+ excited state) of 
^^F is described as a 75 = proton (75 = 51/2 proton) coupled to the ^^O core, which 
has an intrinsic spin I a =0+. The gamma-ray transitions are dominated by the E\ 
multipolarity and by incoming p waves for both states. The Ml and E2 contributions 
amount to less than 0. 1% of the dominant El contribution, as shown in Ref. [11] where 
a potential model was also used. The spin-orbit potential depth —10.0 MeV and the 
spectroscopic factor 1.0 have been used for both states, although 0.9 and 1.0 for the 
ground state and the excited state, respectively, are recommended by Rolfs in Ref. [11]. 
On the other hand, Iliadis et al. [18] recommend values close to unity. The continuum 
states potential depth are set as the same as that of bound states, since no elastic scat- 
tering data is known for this system. Our results are shown in Figure 2 with solid line, 
and as we can see, it reproduces quite well the experimental data. The parameterization 
of the potential model calculations for the ^^0(p,7)^^F capture reactions with the pole 
expression is given by: 

for i^FG5(5/2-h): S{E) 
;t2=0.0139; 

for 17F£5(1/2+): S{E) 

Here also we present the results for the ^lA{n,yf'llgs^\st, ^lA{n,yf'Ligs, ^Be(p,7)^Bgi 

and ^B(p,7)9q„s 

For the ''lA{n,yf''L\gs^ist capture reaction, the gamma-ray transitions are dominated by 
the E\ multipolarity and by incoming s waves and d waves. The Jb = 2+ ground state 
(Jb= l"*" first excited state) of ^Li is described as a Jb = P3/2 neutron interacting with 
the ^Li core, which has an intrinsic spin 4 = 3 /2~. Here we used R() = Rc = Rso = 2.391 



= q^^Ie + 0-3207 + 0.2734 x E, giving S(0)=0.319 and 
= + 2.9640 - 0.1005 x E, with S(0)=10.987 and 
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FIGURE 2. Potential model calculations for the reaction ^Li(n,7)^Li. Experimental data are from Refs. 
[16, 19, 20,21,22]. 



fm for the radius parameters. For the continuum state, the potential depth was obtained 
from the analysis of Nagai [16] which was determined from the experimental scattering 
lengths, see Table-I for the values. The spectroscopic factors adopted where SF{g.s.) = 
0.98 [23] and SF{lst) = 0.48 [16] for the ground and first excited states, respectively. 
The capture to the first excited state contributes to less than 5% of the total cross section. 
The results for this calculation are shown in Figure 2. 

For the ^hi(n,Yf'Ligs capture reaction, the scattering potential parameters for both 
entrance channel spins, s = 5/2+, 3/2+, for the ^Li(2+)+n system were obtained by 
keeping the same volume integral per nucleon, J^/A, as those for the entrance channel 
spins, 5 = 2+, 1+, deduced from the scattering potentials of the ^Li+n system [16]. The 
results are presented in Figure 4. Details of the analysis for this neutron capture reactions 
^Li(n,7)^Lig,5^i^^ and ^Li(n,7)^Ligj are published in Ref. [23]. 

For the ^Be{p,Y)^Bgs capture reaction, the 7g = 2+ ground state of is described as 
a Jb = P3/2 proton coupled to the ^Be core, which has an intrinsic spin /a = 3/2^. For 
this system we adopted a = 0.52 fm and = —9.8 MeV. This is the same set of values 
adopted in Ref. [17]. The gamma-ray transition is dominated by the El multipolarity 
and by incoming s and d waves. Our results are shown in Figure 3. To reproduce the Ml 
resonance we considered Vscat = —38.14 MeV and SF = 0.7 (dashed-dotted line) also 
with the other parameters according to Table I. 

For the ^B(p,7)^Cg.s no data is available. The capture process for this reaction is 
dominated by El transitions from incoming 5-waves to bound p states in the ground 
state of ^C(/e = 3/2^), which is described as a jg = j2 proton coupled to the ^B core, 
which has an intrinsic spin = 2+. The spectroscopic factor has been set to 1.0 as in 
Ref. [30], where several spectroscopic factor values are compared. The results are shown 
in Figure-4. 



FINAL REMARKS 
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FIGURE 3. (Potential model calculations for the reaction ^Be(p,7)^B. The dashed-dotted line is the 
calculation for the Ml resonance at £„„ = 0.63 MeV and the dotted line is for the non-resonant capture. 
Experimental data are from Refs. [24, 25, 26, 6, 27, 28]. The total S factor is shown as a solid line. 




- - U W.^ U.4 U.Q U.O 

E [MeV] Neutron Energy(MeV) 



FIGURE 4. (a) Potential model calculations for the reaction ^B(p,7)^C. The solid circle symbol at £ = 
is from Refs. [29], which were obtained by ANC (extracted from breakup reaction) calculation, (b) Po- 
tential model calculations for the **Li(n,7)'^Li capture reaction. The lower curves labeled (a) correspond to 
the adoption of scattering potential derived from Table-I. Curves labeled (b) correspond to the assumption 
of same potential for the incoming wave function as for the bound state, for s-wave neutron only (dotted 
curve) and s- and d-wave neutrons (solid curve). The uppermost curve (dashed line) corresponds to a 1/v 
fit to the low-energy cross sections. 



In this work, we report the resuhs obtained for the non-resonant part of the neutron 
and proton capture reactions of light nuclei; ^Li(n,7)^Li, ^Li(n,7)^Li, ^Be(p,7)^B and 
^B(p,7)^C in the framework of the potential model. As a test of the model we show the 
resuhs for the i^O(p,7)^^Fg,(5/2+) and ^'^0(p,7)^^Fe^(l/2+) capture reactions. These 
calculations are part of a more extensive series of calculations of neutron and proton 
capture reactions by light radioactive nuclei published elsewhere [31]. 

Although we have obtained the incoming nucleon scattering potentials for the reaction 
of interest from analysis of close systems, or by considering the same as the bound 
state, it would be interesting to obtain such potentials from direct elastic scattering 
measurement as ^Li-i-p and ^B-i-p. A program of investigation for these elastic scattering 



TABLE 1. Wood-Saxon potential parameters used in the calculations. Depths and 
B.E. are in MeV with ro = rso = fcoui = 1-25 fm and a = aso = 0.65 fm, where the 

1/3 

radii are given by /? = roxAj! . The spin-orbit is given by Vjo— 10 MeV. 





B.E. 


Vo (bound) 


SF 


channel spin 


Vo(scat) 


160+p="Fj,,(5/2+) 


0.60 


51.79 


1.0 




51.79 


i6o+p=i^Fi,,(l/2+ 


0.11 


51.42 


1.0 




51.42 




0.14 


41.26 


1.0 




41.26 


^B+p=''q,,(3/2-) 


1.30 


41.97 


1.0 




41.97 


^Li+n=*^Li^,(2+) 


2.03 


46.38 


0.98(15) 


2+,l + 


56.15,46.50 


■^Li+n=^Lii,f(l+) 


1.05 


43.30 


0.48 


2+,l + 


56.15,46.50 


^Li+n=5Li^,(3/2-) 


4.06 


47.82 


0.62(13) 


5/2+, 3/2+ 


58.15,48.15 



experiments at low energy is under way at the Sao Paulo University using the radioactive 
ion beam facility RIBRAS [32]. 
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